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Abstract 

Using 55.8 pb _1 of e + e~ collisions recorded at the ^(3770) resonance with the CLEO-c detector 
at CESR, we determine absolute hadronic branching fractions of charged and neutral D mesons 
using a double tag technique. Among measurements for three D° and six D + modes, we obtain 
reference branching fractions B(D° -> R--k + ) = (3.91 ± 0.08 ± 0.09)% and B{D + -> K-h + tt + ) = 
(9.5 ± 0.2 ± 0.3)%, where the uncertainties are statistical and systematic, respectively. Final state 
radiation is included in these branching fractions by allowing for additional, unobserved, photons 
in the final state. Using a determination of the integrated luminosity, we also extract the cross 
sections a(e + e~ -» D°D°) = (3.60 ± 0.07lg;^) nb and a(e + e- -» D+D~) = (2.79 ± 0.07±g;J2) nb. 
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Absolute measurements of hadronic charm meson branching fractions play a central role 
in the study of the weak interaction because they serve to normalize many D and B meson 
branching fractions, from which elements of the Cabibbo-Kobayashi-Maskawa (CKM) ma- 
trix [lj are determined. For instance, the determination of the CKM matrix element \V c b\ 
from the B — > D*tv decay rate using full D* reconstruction requires knowledge of the D 
meson branching fractions 0. In this Letter, we present charge-averaged branching frac- 
tion measurements of three D° and six D + decay modes: D° — > K~7t + , D° — > K~-k + 7t , 
D° -> K-n + n + n-, D + -> K-n + n + , D + -> R-n+n+n , D + -> f^7r+, L>+ -> ^tt+tt , 
D + — ► iT57r + 7r + 7r _ , and -D + — > K + K~tt + . Two of these modes, _D° — > fC~7r + and 
D + — > fT~7r + 7r + , are particularly important because essentially all other D° and -D 4 " branch- 
ing fractions have been determined from ratios to one of these branching fractions 

To date, the most precise measurements of hadronic D branching fractions are made with 
slow-daughter-pion tagging of D* mesons from Z° decays and from continuum production 
in e + e~ interactions at the T(45) 0, HI- Previously, the MARK III collaboration measured 
hadronic branching fractions at the DD threshold using a double tagging technique which 
relied on fully-reconstructed ^(3770) — > DD decays 0, 0j- This technique obviated the need 
for knowledge of the luminosity or the e + e~ — > DD production cross section. We employ 
a similar technique using CLEO-c data, in a sample roughly six times larger than that of 
MARK III, resulting in precision comparable to the current PDG world averages. 

The data sample we analyze was produced in e + e~ collisions at the Cornell Electron 
Storage Ring (CESR) and collected with the CLEO-c detector. It consists of 55.8 pb -1 
of integrated luminosity collected on the ^(3770) resonance, at a center-of-mass energy 
-^cm — 3773 MeV. At this energy, no additional hadrons accompanying the DD pairs are 
produced. Reconstruction of one D or D meson (called single tag or ST) tags the event as 
either D°D° or D + D~ . For a given decay mode i, we measure independently the D and D 
ST yields, denoted by N{ and N{. We determine the corresponding efficiencies, denoted by 
€i and 6j, from Monte Carlo simulations. Thus, Ni = eiBiN D p and Ni = €iBiN D £), where 
Ei is the branching fraction for mode i, assuming no CP violation, and N d q is the number 
of produced DD pairs. Double tag (DT) events are the subset of ST events where both 
the D and D are reconstructed. The DT yield for D mode i and D mode j, denoted by 
Ny, is given by Nij = €ijBiBjN D Q, where is the DT efficiency. As with ST yields, the 
charge conjugate DT yields and efficiencies, Nji and e 3 -$, are determined separately. Charge 
conjugate particles are implied, unless referring to ST and DT yields. 

The Bi can be determined from the DT yield N^ and the corresponding ST yield Nj 
via Bi = [Nij/Nj] x [ej/e^]. Similarly, we have N D £, = [NiNj/Nij] x [cij/(€i€j)]. Because 
tij 6i€j, the branching fractions thus obtained are nearly independent of the tag mode 
efficiency, and N D d is nearly independent of all efficiencies. We extract branching fractions 
and N D £) by combining ST and DT yields with a least squares technique. Although the 
D° and D + yields are statistically independent, systematic effects and misreconstruction 
resulting in crossfeed introduce correlations among their uncertainties. Therefore, we fit D° 
and D + parameters simultaneously, including in the x 2 statistical and systematic uncer- 
tainties and their correlations for all experimental inputs |8|. Thus, yields, efficiencies, and 
backgrounds are treated uniformly, and the statistical uncertainties on Bi and N D ^ include 
the correlations among Ni, Nj, and N^. Also, in the above efficiency ratios most systematic 
uncertainties are correlated between ST and DT efficiencies, so their effects largely cancel. 

The CLEO-c detector is a modification of the CLEO III detector j^, 0, u^, in which 
the silicon-strip vertex detector was replaced with a six-layer vertex drift chamber, whose 
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wires are all at small stereo angles to the beam axis lil . The charged particle trac king 
system, consisting of this vertex drift chamber and a 47-layer central drift chamber [lOj 
operates in a 1.0 T magnetic field, oriented along the beam axis. The momentum resolution 
of the tracking system is approximately 0.6% at p = 1 GeV/c. Photons are detected in an 
electromagnetic calorimeter, composed of 7800 CsI(Tl) crystals 0, which attains a photon 
energy resolution of 2.2% at E 1 = 1 GeV and 5% at 100 MeV. The solid angle coverage for 
charged and neutral particles of the CLEO-c detector is 93% of 4tt. We utilize two particle 
identification (PID) devices to separate from t[ ± : the central drift chamber, which pro- 
vides measurements of ionization energy loss (dE/dx), and, surrounding this drift chamber, 
a cylindrical ring- imaging Cherenkov (RICH) detector [111 ], whose active solid angle is 80% 
of An. The combined dE/dx-KlCB. PID system has a pion or kaon efficiency > 90% and 
a probability of pions faking kaons (or vice versa) < 5%. The response of the experimen- 
tal apparatus is studied with a detailed GEANT-based [l3| Monte Carlo simulation of the 
CLEO detector for particle trajectories generated by EvtGen and final state radiation 
(FSR) predicted by PHOTOS [15]. Simulated events are processed in a fashion similar to 
data. The data sample's integrated luminosity (£) is measured using e + e~ Bhabha events 
in the calorimeter [16], where the event count normalization is provided by the detector 
simulation. 

Charged tracks are required to be well-measured and to satisfy criteria based on the track 
fit quality. They must also be consistent with coming from the interaction point in three 
dimensions. Pions and kaons are identified by consistency with the expected dE/dx and 
RICH information, when available. We form tt° candidates from photon pairs with invariant 
mass within 3 standard deviations (a), with a rs 5-7 MeV/c 2 depending on photon energy 
and location, of the known 7r° mass. These candidates are then fit kinematically with their 
masses constrained to the known tt° mass. The Kg candidates are selected from pairs of 
oppositely-charged and vertex-constrained tracks having invariant mass within 12 MeV/c 2 , 
or roughly 4.5a, of the known K s mass. 

We identify D meson candidates by their invariant masses and total energies. We calculate 
a beam-constrained mass by substituting the beam energy, Eq, for the measured D candidate 



energy: Mc 2 = \J E 2 — p^c 2 , where is the D candidate momentum. Performing this 
substitution improves the resolution of M by one order of magnitude, to about 2 MeV/c 2 , 
which is dominated by the beam energy spread. We define AE = Ed — E , where Ed is 
the sum of the D candidate daughter energies. For final states consisting entirely of tracks, 
the AE resolution is 7-10 MeV. A 7r° in the final state degrades this resolution by roughly 
a factor of two. We accept D candidates with M greater than 1.83 GeV/c 2 and with mode- 
dependent AE requirements of approximately 3a. For both ST and DT modes, we accept 
at most one candidate per mode per event. In ST modes, the candidate with the smallest 
AE is chosen, while in DT modes, we take the candidate whose average of D and D M 
values, denoted by M, is closest to the known D mass. 

We extract ST and DT yields from M distributions in the samples described above. We 
perform unbinned maximum likelihood fits in one and two dimensions for ST and DT modes, 
respectively, to a signal shape and one or more background components. The signal shape 
includes the effects of beam energy smearing, initial state radiation, the line shape of the 
^(3770), and reconstruction resolution. The background in ST modes is described by an AR- 
GUS function 0], which models combinatorial contributions. In DT modes, backgrounds 
can be uncorrelated, where either the D or D is misreconstructed, or correlated, where all the 
final state particles in the event are correctly reconstructed but are mispartitioned among 
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TABLE I: Single tag data yields and efficiencies and their statistical uncertainties. 



D or 


D Mode 




Yield (10 3 ) Efficiency (%) 


r-iO 

L> v — 


r ' _|_ 






5.11 ± 0.07 


64.6 ± 0.3 


w — 


> if + 7T 






5.15 ± 0.07 


65.6 ± 0.3 


7-1O 

D — 


> if 7T + 


7T° 




9.51 ± 0.11 


31.4 ± 0.1 


JJ — 




7T° 




9.47 ± 0.11 


31.8 ± 0.1 


7-1O 


> K 7T + 


^+ 


7T 


7.44 ± 0.09 


43.6 ± 0.2 


iJ u — 


tv-4- — 
> if + 7T 


7T~ 


^+ 


7.43 ± 0.09 


43.9 ± 0.2 


D + - 


-> if 7T H 


-7T+ 


7.56 ± 0.09 


50.7 ± 0.2 


D - 


7- -J- 1 


~7T 




7.56 ± 0.09 


51.3 ± 0.2 


£> + - 






-7T° 


2.45 ± 0.07 


25.7 ± 0.2 


D~ - 


-» K + it- 




'7T° 


2.39 ±0.07 


25.7 ±0.2 


D+ - 








1.10 ±0.04 


45.5 ±0.4 


D~ - 


^K° S TT- 






1.13 ±0.04 


45.9 ±0.4 


D+ - 


^if°7T+ 


7T° 




2.59 ±0.07 


22.4 ±0.2 


D~ - 


-> KgTT~ 


7T° 




2.50 ±0.07 


22.4 ±0.2 


D+ - 


+ K° S *+ 


7T+ 


7T~ 


1.63 ±0.06 


31.1 ±0.2 


D~ - 




7T~ 


7T+ 


1.58 ±0.06 


31.3 ±0.2 


D+ - 


■* K + K 


~7T 


+ 


0.64 ±0.03 


41.4 ±0.5 


D~ - 


■* K+K 


~7T 




0.61 ±0.03 


40.8 ±0.5 



the D and D. In fitting the two-dimensional M(D) versus M(D) distribution, we model the 
uncorrelated background by a pair of functions, where one dimension is an ARGUS func- 
tion and the other is the signal shape. We model the correlated background by an ARGUS 
function in M and a Gaussian in the orthogonal variable, which is [M(D) — M(D)}/2. 

Table |U gives the 18 ST data yields and efficiencies determined from simulated events. 
Figure ^ shows the M distributions for the nine decay modes with D and D candidates 
combined. Overlaid are the fitted signal and background components. We also measure 
45 DT yields in data and determine the corresponding efficiencies from simulated events. 
Figure 121 shows M(D) for all modes combined, separated by charge. We find total DT yields 
of 2484 ±51 for D° and 1650 ±42 for D+. Because of the cleanliness of the DT modes, their 
statistical yield uncertainties are close to yNij- 

Using a missing mass technique, we measure efficiencies for reconstructing tracks, Kg 
decays, and n° decays in both data and simulated events. We fully reconstruct ?/>(3770) — > 
DD, ip{2S) — > J/ipir + iT~, and ij){2S) — > J/ip7i°7i° events, leaving out one particle, for which 
we wish to determine the efficiency. The missing mass of this combination peaks at the 
mass of the omitted particle, whether or not it is detected. Then, the desired efficiency 
is the fraction of this peak with this particle correctly reconstructed. For tracks and K s 
candidates, we find good agreement between efficiencies in data and simulated events. For 7r° 
candidates, we correct the simulated efficiencies by 3.9%, which is the level of disagreement 
with data found in this study. The relative uncertainties in these determinations, 0.7% 
per track, 3.0% per K s , and 2.0% per n°, are the largest contributions to the systematic 
uncertainties. 
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FIG. 1: Semilogarithmic plots of ST yields and fits, with D and D combined for each mode. Data 
are shown as points with error bars. The solid lines show the total fits and the dashed lines the 
background shapes. The high mass tails on the signal are due to initial state radiation. 



We study the simulation of the PID efficiencies using decays with unambiguous parti- 
cle content, such as D° — > K^n + 'K~ and D + — > K~tt + tt + . We find a need to correct the 
simulated efficiencies by 0.3% for ir^ and 1.3% for K , and we apply associated system- 
atic uncertainties of the same size. Other sources of efficiency uncertainty include: the AE 
requirements (1.0-2.5%), for which we examine AE sidebands; modeling of particle multi- 
plicity and detector noise (0.2-1.3%); and modeling of resonant substructure in multi-body 
modes (0.4-1.5%), which we assess by comparing simulated momentum spectra to those in 
data. We also include additive uncertainties of 0.5% to account for variations of yields with 
fit function. Smaller systematic uncertainties arise from online and offline filtering (0.4%), 
simulation of FSR (0.5% per D or D), and the assumed width of the ^(3770) in the M 
signal shape (0.6%). The effect of quantum correlations between the D° and D° states 
appears through D°-D° mixing and through doubly Cabibbo-suppressed decays [l8|. The 
former contribution is limited by available measurements to be less than O(10~ 3 ) and is, 
therefore, neglected in this analysis. The latter contribution is addressed with a systematic 
uncertainty due to the unknown phase of interference in neutral DT modes between the 
Cabibbo-favored amplitude and the amplitude for doubly Cabibbo-suppressed transitions in 
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FIG. 2: DT yields and fits projected onto the D° and D + axes and summed over all modes. Data 
are shown as points with error bars. The solid lines show the total fit and the dashed lines the 
background shapes. 

both D° and D° (0.8%). 

The branching fraction fitter 0] takes these systematic uncertainties as input, along with 
ST and DT yields and efficiencies, crossfeed probabilities among the modes, background 
branching fractions and efficiencies, and statistical uncertainties on all of these measure- 
ments. The estimated crossfeed and background contributions produce yield adjustments 
of (9(1%). Their dependence on the fit parameters is taken into account both in the yield 
subtraction and in the \ 2 minimization. We validated the analysis technique, including the 
branching fraction fit, by studying simulated DD events in a sample 50 times the size of our 
data sample. We reproduced the input branching fractions with biases due to our procedures 
that were less than one-third of the statistical errors on the data and consistent with zero. 

The results of the data fit are shown in Table ITT1 The x 2 °f the fit is 28.1 for 52 degrees 
of freedom, corresponding to a confidence level of 99.7%. To obtain the separate contribu- 
tions from statistical and systematic uncertainties, we repeat the fit without any systematic 
inputs and take the quadrature difference of uncertainties. All nine branching fractions are 
consistent with, and most are higher than, the current PDG averages 0]. In the D candi- 
date reconstruction, we do not explicitly search for FSR photons. However, because FSR is 
simulated in the samples used to calculate efficiencies, our measurements represent inclusive 
branching fractions for signal processes with any number of photons radiated from the final 
state particles. If no FSR were included in the simulations, then all the branching fractions 
would change by Afsr in Table |H] 

The correlation coefficient, p, between N D o£,o and N D + D - is 0.07, and each is essentially 
uncorrelated with branching fractions of the other charge. Correlations among branching 
fractions are in the range 0.2-0.7. In the absence of systematic uncertainties, there would 
be almost no correlation between the charged and neutral D parameters. 
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TABLE II: Fitted branching fractions and DD pairyields, along with the fractional FSR correc- 
tions and comparisons to the Particle Data Group [3j fit results. Uncertainties are statistical and 
systematic, respectively. 



D Decay Mode Fitted B (%) 


PDG B (%) 


Afsr 


K--K+ 


3.91 ±0.08 ±0.09 


3.80 ± 0.09 


-2.0 


% 


K-TT+TT 


14.9 ±0.3 ±0.5 


13.0 ±0.8 


-0.8 


% 


K~TT + ir + TT~ 


8.3 ±0.2 ±0.3 


7.46 ± 0.31 


-1.7% 




9.5 ±0.2 ±0.3 


9.2 ±0.6 


-2.2 


% 


K-ir+ir+ir 


6.0 ±0.2 ±0.2 


6.5 ± 1.1 


-0.6 


% 


K° S 7T + 


1.55 ±0.05 ±0.06 


1.41 ±0.10 


-1.8 


% 




7.2 ±0.2 ±0.4 


4.9 ±1.5 


-0.8 


% 




3.2 ±0.1 ±0.2 


3.6 ±0.5 


-1.4 


% 


K+K-1T+ 


0.97 ±0.04 ±0.04 


0.89 ±0.08 


-0.9 


% 


DD Yield 


Fitted Value 


Afsr 


N D o£)0 
N D+D - 


(2.01 ±0.04 ±0.02) x 10 5 
(1.56 ±0.04 ±0.01) x 10 5 


-0.2 
-0.2 


% 
% 



TABLE III: Fitted ratios of branching fractions to the reference branching fractions TZo = B(D° — > 
K~tt + ) and 1Z+ = B{D + — > K~7r + ir + ), along with the fractional FSR corrections and comparisions 
to the Particle Data Group jjj fit results. Uncertainties are statistical and systematic, respectively. 



D Decay Mode 


Fitted B/K 0/+ 


PDG B/K Q/+ 


Afsr 


K-7T+TT 


3.82 ±0.05 ±0.10 


3.42 ±0.22 


+1.2% 




2.12 ±0.03 ±0.06 


1.96 ±0.08 


±0.3% 


K-ir+ir+ir 


0.634 ±0.014 ±0.018 


0.70 ±0.12 


±1.7% 


K° S TT+ 


0.162 ±0.004 ±0.006 


0.153 ±0.003 


±0.4% 


K° S TT + TT° 


0.753 ±0.016 ±0.039 




±1.4% 




0.336 ± 0.009 ± 0.014 


0.39 ± 0.05 


±0.8% 


K+K--K+ 


0.102 ±0.004 ±0.003 


0.097 ± 0.006 


±1.3% 



We also compute ratios of branching fractions to the reference branching fractions, shown 
in Table ILTTl These ratios have higher precision than the individual branching fractions, and 
they also agree with the PDG averages. Without FSR corrections to the efficiencies, all 
seven ratios would be 0.3% to 1.7% higher. 

We obtain the e + e~ — > DD cross sections by scaling N D of)0 and N^+d- by the lumi- 
nosity, which we determine to be C = (55.8 ± 0.6) pb _1 . Thus, at E cm = 3773 MeV, 
we find peak cross sections of a{e + e~ -> D°D°) = (3.60 ± 0.07±g;^) nb, a{e + e~ -> 
D+D-) = (2.79 ± 0.07±8;J2) nb, a{e + e~ -> DD) = (6.39 ± O.lOlg;^) nb, and cr(e+e" -> 
D+D )/a(e + e- -> D°D°) = 0.776 ± 0.024ig;gJ|, where the uncertainties are statistical and 
systematic, respectively. In addition to the systematic uncertainties on N d oqo, N d + d -, and 
the luminosity, the cross section systematic uncertainties also include the effect of E cm vari- 
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ations with respect to the peak. We account for the correlation between the charged and 
neutral cross sections in computing the uncertainty on the total cross section. Our measured 
cross sections are in good agreement with BES 0] and higher than those of MARK III [?| . 

In summary, we report measurements of three D° and six D + branching fractions and the 
production cross sections a(D°D°), a(D + D~), and a(DD) using a sample of 55.8 pb _1 of 
e + e~ — > DD data obtained at E cm = 3773 MeV. We find branching fractions in agreement 
with, but somewhat higher, than those in the PDG summary. We note that, unlike the 
branching fractions used in the PDG averages, our measurements are corrected for FSR. 
Not doing so would lower our branching fractions by 0.6% to 2.2%. With our current data 
sample, the statistical and systematic uncertainties are of comparable size. Many of the 
systematic uncertainties, such as those for tracking and particle identification efficiencies, 
will be improved with larger data samples. 

We gratefully acknowledge the effort of the CESR staff in providing us with excellent 
luminosity and running conditions. This work was supported by the National Science Foun- 
dation and the U.S. Department of Energy. 
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